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Abstract. Particle acceleration at supernova remnant (SNR) shock waves is regarded as the most probable mechanism for 
providing Galactic cosmic rays at energies below 10^^ eV. The Galactic cosmic ray hadron component would in this picture 
result from the injection of relativistic particles from many SNRs. It is well known that the superposition of individual power 
law source spectra with dispersion in the spectral index value, which behaviour is observed in the synchrotron radio spectra of 
shell SNR, displays a positive curvature in the total spectrum and in particular shows a hardening at higher energies. Recent ob- 
servations made with the EGRET instrument on the Compton Gamma-Ray Observatory of the diffuse Galactic 7-ray emission 
reveal a spectrum which is incompatible with the assumption that the cosmic ray spectra measured locally hold throughout the 
Galaxy: the spectrum above 1 GeV, where the emission is supposedly dominated by vr^-decay, is harder than that derived from 
the local cosmic ray proton spectrum. We demonstrate that in case of a SNR origin of cosmic ray nucleons part of this 7-ray 
excess may be attributed to the dispersion of the spectral indices in these objects. In global averages, as are 7-ray line-of-sight 
integrals, this dispersion leads to a positive curvature in the composite spectrum, and hence to modified 7r"-decay 7-ray spectra. 
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1. Introduction 

Particle acceleration at supernova remnant shock waves is re- 
garded as the most probable mechanism for providing Galactic 
cosmic rays at energies below 10^'' eV (for a review see 
Blandford & Eichler 1987). The recent detections of non- 



thermal X-ray synchrotron radiation from the four supernova 
remnants SN1006 (Koyama et al.|l995|), RX J 17 13.7-3946 



(Koyama et al. 1997), Slane et al. 1999), Cas A (Allen et al 



1997), and RCW86 (Borkowski et al. 2001) and the subsequent 



detections of SN1006 (Tanimori et al. 1998) and RX J1713.7 



3946 (Muraishi et al. 200C ) at TeV energies support the hypoth- 



esis that at least Galactic cosmic ray electrons are accelerated 
predominantly in SNR. To date, there is still no unambiguous 
proof that cosmic ray nucleons are similarly produced in SNR. 

Whatever the nature of their sources, upon escape from 
there the cosmic ray nucleons would diffusively propagate 
through the interstellar medium, where they can be either di- 
rectly measured or indirectly traced by means of 7-ray ob- 
servations. Recent observations made with the EGRET instru- 
ment on the Compton Gamma-Ray Observatory of the diffuse 
Galactic 7-ray emission reveal a spectrum which is incompat- 
ible with the assumption that the cosmic ray spectra measured 



locally hold throughout the Galaxy (Hunter et al. 1997). The 



spectrum observed with EGRET below 1 GeV is in accord 
with, and supports, the assumption that the cosmic ray spec- 
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tra and the electron-to-proton ratio observed locally are uni- 
form, however, the spectrum above 1 GeV, where the emission 
is supposedly dominated by 7r° -decay, is harder than that de- 
rived from the local cosmic ray proton spectrum. 

In a recent paper Pohl & Esposito ( 1998 ) demonstrated that 
if the sources of cosmic rays are discrete, as are SNR, then 
the spectra of cosmic ray electrons would vary and the locally 
measured electron spectrum would not be representative for the 
electron spectra elsewhere in the Galaxy, which could be sub- 
stantially harder than the local one. These authors have shown 
that the observed excess of 7-ray emission above 1 GeV can in 
fact be explained as a correspondingly hard inverse Compton 
component, provided the bulk of cosmic ray electrons is pro- 
duced in SNR. 

In the following we will demonstrate that in case of a SNR 
origin of cosmic ray nucleons part of the 7-ray excess may also 
be attributed to the dispersion of the spectral indices in these 
objects. In global averages, as are 7-ray line-of-sight integrals, 
this dispersion leads to a positive curvature in the composite 
spectrum, and hence to modified tt" -decay 7-ray spectra. 

2. The dispersion of cosmic ray spectra in SNR 

The synchrotron radio spectra from shell S NR in dicate power 
law sp ectra oc v'" (Clark & Caswell |l976l Milne |l 979 , 
Green 2001 ) where the spectral index distribution from the 
sample of shell supernova remnants (SNR) has a mean value 
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of about < a >~ 0.5 and a significant standard deviation 
a a- From synchrotron radiation theory this imphes a power 
law energy distribution of the radiating relativistic electrons 
with mean spectral index Se = 1 + 2 < a >~ 2.0 and 
dispersion CTe — 2a a- Because the age of the SNRs is much 
smaller than the characteristic radiative loss times of both cos- 
mic ray nucleons and GeV electrons in the remnant, and be- 
cause both the particle acceleration processes and the spatial 
propagation scale with rigidity, the relativistic hadrons should 
have the same momentum spectrum as the relativistic electrons. 
After leaving their source these particles propagate in the inter- 
stellar medium by momentum-dependent spatial diffusion with 
a diffusion coefficient k oc p^, with b ~ 0.6 as inferred from 
the measurement of secondary to primary cosmic ray elements. 
Consequently, the source spectrum of cosmic ray hadrons 



mc 



(1) 



is a power law in momentum per nucleon with a spectral index 
having a mean value < s >= 1 + 6 + 2 < a >~ 2.7 and 
dispersion a — Ue- 

The Galactic cosmic ray hadron component results from the 
injection of relativistic particles from many SNRs, especially 
in the inner part of the Galaxy. In discussing the pros and cons 
of the Galactic origin of cosmic rays, Brecher and Burbidge 
(1972) noted that the superposition of individual power law 
source spectra with dispersion in the spectral index value dis- 
plays a positive curvature in the total cosmic ray spectrum and 
in particular shows a flattening at higher energies, simply be- 
cause those sources with the smallest spectral index dominate 
the total spectrum at large energies. If we represent the distri- 
bution of hadron spectral indices by the Gaussian 



n{s) 



1 
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exp 
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(2) 



we obtain from Eq. (|l]) for the averaged hadron spectrum 

<N{p)>= / dsN{p,s)n{s) 



\mcJ 



-<s> + V In 



(3) 



By using N{p)dp — N{'-f)dj and p = mcy^j'^ — 1 we can 
determine from Eq. the corresponding differential hadron 
spectrum with respect to the hadron Lorentz factor as 

N{-f) = Name 7(7^ - 1)^^^+^ 1.1(7"-!) 

With increasing hadron Lorentz factor 7 the energy spectrum 
(jsj) flattens, i.e. it hardens. Calculating the slope yields 

_ _ In < iV(7) > 



1117 

< s > +1 ln(72 - 1) 
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At relativistic energies Eq. (||) approaches 
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3. Pion decay gamma rays 

For the tt" 27 decay the omnidirectional (i.e. integrated 
over the whole solid angle) differential 7-ray source function 
g^o [E-f ; r) at the position r = {l,b,r) in space is related to the 
omnidirectional differential neutral pion source function 



with the lower integration boundary 



C/(i?^) 
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in terms of the dimensionless 7-ray photon energy 

E-y 



0.135 GeV 



(7) 



(8) 



(9) 



in units of the pion rest energy. 

The differential photon number flux of 7r°-decay 7-rays 
from a direction (1, b) is given by the line of sight integral of 
the source function (ITh 



dN^{E^;l,b) 
dt dE-y dQ 



1 
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dr q^o {E^;r) 
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3. 1. Pion source spectrum 

The pion source spectrum has been frequently calculated on 
the basis of the inclusive cross section and a delta-functional 



for the pion spectrum (e.g Mannheim & Schlickeiser |1994[ ). As 
we show in the Appendix this approximation may be useful at 
higher energies, but it is certainly inappropriate at energies not 
very far from the threshold. 

Therefore in this section the pion source spectrum required 
in Eq. (0) is calculated using cosmic hadron distribution func- 
tion (|) as input f or the Monte-Car lo cod e DTUNUC (V2.2) 
(Mohring & Ranft|l99l| , Ranft et al. |l994i Ferrari et al. |l996a , 
Engel et al. 1997), which is based on a dual parton model 



(Capella et al. 1994 ). This MC model for hadron-nucleus and 
nucleus-nucleus interactions includes various modern aspects 
of high-energy physics and has been successfully applied to 
the description of hadron production in high-energy collisions 
(Ferrari et al. |l996b|, Ra nft & Roesler |1994 Mohring et al. 



1993,Roesleretal. 1998) 



The 7-ray spectra thus derived can be compared with the 
observed EGRET spectra from the inner Galaxy. Here we use 
the data of all viewing periods of phases 1-4, corresponding 
to observations between 1991 April and 1995 October From 
the observed intensity we subtract the extragalactic background 
(Sreekumar et al. 1998) and the point-spread functions of all 
sources in the Third EGRET Catalogue (Hartman et al. 1999). 



The intensity spectrum has additional associated uncertainties, 
which arise from the subtraction of the extragalactic back- 
ground and the sources. A systematical error of 10% is assumed 
to account for uncertai nties in the spark chamber efficiency cor- 
rection (Esposito et al. 1999). 

The pion decay spectra are compared with the intensity of 
the diffuse Galactic 7-ray emission from the inner Galaxy in 
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Fig. 1. The observed intensity spectrum from the inner Galaxy 
shown in comparison with 7r° -decay spectra with (bottom 
panel) and without (top panel) dispersion in the cosmic ray 
spectrum. The solid lines display the 7r°-decay spectra, and the 
dashed lines are the total 7-ray spectra including the leptonic 
contribution, which here is simply given as a power law oc i?^^ 
with an intensity determined by the data for E < 100 MeV. 
Whereas in the top panel, i.e. without dispersion, the GeV ex- 
cess can be clearly seen, the bottom panel proves that a disper- 
sion with <T = 0.25 can in fact explain the data and remove the 
GeV excess. 

Fig. |l[ There the leptonic contribution to the diffuse 7-ray in- 
tensity has been approximated by a simple power law c>c E^'^ 
with an intensity E^Ii = 0.01 McV/cm'^/sec/sr, adjusted 
to reproduce the data for E < 100 MeV. It can be seen that 
a dispersion in the index of the cosmic ray source spectra of 
a ~ 0.25 is sufficient to explain the observed intensity spec- 
trum and remove the GeV excess. We can therefore conclude 
that if dispersion is to be made responsible for the GeV excess, 
it must be at a level of ct ~ 0.25. 

4. Compatibility withi upper iimits at liiglier 
gamma ray energies 

A dispersion in the cosmic ray source spectra - and thus in the 
cosmic ray flux throughout the Galaxy - has impact also on the 
resultant 7-ray spectra at higher energies, both for the diffuse 
emission and for individual supernova remnants. In this section 
we will investigate whether a spectral dispersion of cr = 0.25 



is compatible with existing upper limits for 7-ray emission in 
the TeV energy range. 

To date, observations of TeV 7-ray emission from indi- 
vidual SNR have yielded only a few detections, one of which 
(Cas A, Aharonian et al. 2001a) is presumably (Atoyan et al. 
2000) and the others are clearly caused by leptonic emission. 
Although deep surveys have been performed, no hadronic TeV 
7-ray emission has been unambiguously detected from SNR 
(Buckley et al. 1998). These observational results have been a 



severe constraint for simple and conventional models of cos- 
mic ray acceleration in SNR. If there were a dispersion in the 
cosmic ray spectra in SNR, some, i.e. those with a hard spec- 
trum, would be expected to be very prominent TeV sources. 
Therefore the observational constraints are even more severe if 
there is a noteworthy dispersion in the particle spectra. 

However, these constraints rely on the assumption that the 
power law spectra of cosmic rays in the SNR persist to very 
high energies like lO^''^ — 10^^ eV. This may not be a valid 
assumption; in fact more detailed models of particle accel- 
eration including non-linear effects indicate that the resultant 
particle spectra may show cut-offs at energies much smaller 
than 10^'' eV (Baring et al. |1999| ). There is also observational 



evidence that SNR do not produce straight power laws up to 
10^^ eV. For all SNR the X-ray flux, synchrotron or not, is less 
than the extrapolated radio synchrotron spectrum. Since many 
of the sources, in particular the five historical remnants, are too 
young for the electron spectra to be hmited by energy losses, 
acceleration cut-offs must occur at electron energies of 100 
TeV or less (Reynolds & Keohane 1999). At these energies the 
acceleration process should operate similarly for electrons and 
ions, and similar cut-offs must be expected for the cosmic ray 
nuclei, which would be intrinsic to the actual acceleration pro- 
cess. Once a remnant is in the Sedov phase, the maximum en- 
ergy of accelerated particles rises very little. Cas A is probably 
still in the free expansion phase, but proper-motion measure- 
ments of Kepler, Tycho and SN1006 indicate that these objects 
are well in the Sedov dynamics (Moffett, Goss and Reynolds 



1993 ). These remnants are therefore unlikely to accelerate cos- 
mic rays to energies higher than 10^** eV at any time. 

We have made a statistical argument in this paper and there- 
fore may not discuss individual sources, but generally speaking 
it would seem that a dispersion in the spectra of cosmic rays in 
SNR is not in conflict with the existing upper limits for TeV 7- 
ray emission, when an upper cut-off in the cosmic ray spectra is 
taken into account. The same argument applies to the Galactic 
diffuse 7-ray emission at TeV energies as we discuss in some- 
what more detail. 

The Whipple collaboration has recently published an upper 
limit for 7-ray emission at 500 GeV from the Galactic plane, 
based on observations in a small area in the plane defined by 
38.5° < I < 41.5° and -2° < & < 2° (LeBohec et al. pOOOl). 



An upper limit for 7-ray emission at 1 TeV from the Galactic 
plane (38° < I < 43°, -2° < 6 < 2°) h as been reported by the 
HEGRA collaboration (Aharonian et al. 2001b| ). Together with 
the 7-ray intensity in the GeV range these observations con- 
strain models of the GeV excess. Composite spectra of diffuse 
Galactic 7-ray emission from 1 GeV to 500 GeV can be com- 
piled, however, they are subject to systematic uncertainties for 
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a number of reasons. The numerical value of this upper limit at 
500 GeV depends to some extent on the 7-ray spectrum itself. 
Also the observed portion of sky is so small that the finite extent 
of the EGRET point-spread function becomes non-negligible. 
To derive such a composite spectrum of diffuse 7-ray emission 
we will therefore use EGRET data for energies above 1 GeV 
only, for which the 68% containment radius is < 1.7°. To ac- 
count for the effect of the point-spread function and to increase 
the number of counts, we will also derive the EGRET inten- 
sity spectrum from an area slightly larger than that observed 
with Whipple, that is 37° < / < 43° and -2.5° <b< 2.5°. 
When comparing the 7r°-decay spectrum with the composite 
spectrum of diffuse Galactic 7-ray emission, it may be suf- 
ficient to use an analytical approximation for the differential 
pion production cross section, because the uncertainty thus im- 
posed is probably not larger than the systematic uncertainty in 
the data, but the procedure is substantially simplified. Our an- 
alytical approximation for the calculation of the pion source 
spectra is described in Appendix A. 

Given the 7-ray intensity at GeV energies, the model inten- 
sity at 500 GeV depends on a) the mean proton spectral index 
(fixed to < s >— 2.7 in this study), b) a possible dispersion 
in the spectral indices of the cosmic ray sources {a = 0.25 
is required to account for the GeV excess), and c) the limit- 
ing energy to which the power law spectra of the cosmic ray 
sources persist. Since we want to investigate whether or not the 
required dispersion cr = 0.25 is compatible with the Whipple 
upper limit, the only number which we may vary is the high 
energy cut-off in the proton spectrum. Therefore the question 
is: what is the high energy cut-off 7max required to satisfy the 
Whipple upper limit at 500 GeV, given the observed GeV in- 
tensity, < s >= 2.7, and a = 0.25. Then we may ask whether 
or not 7inax thus derived is reasonable in view of the cosmic ray 
all-particle spectrum and the observed composition near and at 
the knee. 

As can be seen in Fig. ^ a dispersion of ct = 0.25 clearly 
violates the upper limit at 500 GeV derived by the Whipple 
team, if 7,„ax = 10^, but the upper limit is clearly satisfied, 
if 7max = 10'^. Therefore, if a spectral dispersion exists at the 
level required to explain the GeV excess, there must be a high 
energy cut-off in the cosmic ray source spectra somewhere be- 
tween 7max = 10^ and 7inax = 10^. 

5. Discussion 

In this paper we have calculated the high energy diffuse 
Galactic 7-ray emission produced by the hadronic component 
of cosmic rays under the assumption that the interstellar cosmic 
ray spectrum is a superposition of individual power-law spec- 
tra with a dispersion in the spectral index as a result of cosmic 
ray production in SNR. We have shown that a diffuse 7-ray 
spectrum thus derived can in fact explain the GeV excess, pro- 
vided the dispersion in the individual SNR production spectral 
indices is cr = 0.25. 

We have then investigated under what assumptions a com- 
posite cosmic ray spectrum with cr = 0.25 would be compati- 
ble with data in the TeV energy range both for individual SNR 
and for the diffuse emission. It is found that the high energy 
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Fig. 2. Composite spectrum between 1 GeV and 1 TeV of 
the diffuse Galactic 7-ray emission from the Galactic plane at 
I sa 4 0° base d on EGRET, Whipple and HEGRA (Aharonian 
et al. ^001b[ ) data. The upper limit at 500 GeV (Whipple), 
which is at the 99.9% confidence level, depends on the 7-ray 
spectrum. The lower mark applies for a 7-ray spectral index 
s — 2.0 whereas the upper mark is appropriate for a spec- 
tral index s = 2.6. The limit at 1 TeV (HEGRA), which is 
a the 99% confidence level, is plotted for two different ways 



of background subtraction (see Aharonian et al. 2001b). The 
data are compared with Tr'^-decay spectra calculated for a cos- 
mic ray spectrum produced by individual sources with mean 
spectral index < s >= 2.7 and dispersion cr = 0.25. The 
solid line refers to a high energy cut-off in the cosmic ray pro- 
ton spectrum at 7max = 10^. It should be compared with the 
lower mark for the upper limit at 500 GeV which it clearly 
and significantly exceeds. It also exceeds the HEGRA upper 
limits at 1 TeV for both ways of background subtraction. For 
7max = 10^ the dashed line shows no contradiction with the 
Whipple and HEGRA data, when compared to the upper mark 
for the upper limit at 500 GeV. 



data would require the existence of a high energy cut-off in 
the cosmic ray source spectra at proton energies not higher 
than somewhere between Lorentz factors 7inax — 10^ and 

7max = 10^. 

The hypothesis presented here provides a possible explana- 
tion for the GeV excess, as do other models. It has a number 
of consequences for observable quantities which can be used to 
test the viability of the scenario and which we want to discuss 
in this section. 

Is a dispersion of ct = 0.25 compatible with the radio 
synchrotron spectra of SNR? The following discussion appUes 
only to shell-type SNR which in fact have an average syn- 
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chrotron power index a ~ 0.5 (Green 2001). This translates 
to an electron spectral index s ~ 2. For various reasons the un- 
certainties in the derived radio spectral indices of SNR are of- 
ten large, but nevertheless some shell-type SNR show spectral 
indices which are significantly different from 0.5, either flatter 
or steeper than the average spectrum. The actual distribution of 
observed spectral indices results from both an intrinsic disper- 
sion and the observational uncertainties, and it is not a trivial 
task to deduce the intrinsic dispersion. We have taken the ra- 
dio spectral index data for Galactic shell SNR given in Fig. 5 of 
Green (2001 1 and have performed a x^-test to see whether the 
sample of measurements is consistent with a uniform source 
spectral index. The answer is clearly "no", the reduced is 
around 12. We have then approximated the effect of a distri- 
bution of source spectral indices by adding a dispersion CTq, to 
the measurement errors. A reduced x^ of unity results when 
the dispersion is Oa — 0.10, then also < a >= 0.53 ± 0.02. 
Taking out Cas A, which has a soft spectrum with very small 
error bars, gives <Ja = 0.085 and < a >~ 0.52 ± 0.02. This 
analysis indicates that a dispersion in the radio spectral indices 
of shell SNR exists, but only at a level of CTq = 0.085 - 0.1, 
corresponding to ag = 0.17 — 0.2 in the cosmic ray electron 
spectra. This falls only marginally short of the a — 0.25 which 
is required to explain the GeV excess. Therefore, given the un- 
certainties, we conclude that a = 0.25 is compatible with the 
radio synchrotron spectra of SNR. 

Why don't we observe a dispersion in the local cosmic ray 
spectrum? Here it is useful to calculate the number of SNR 
which would actually contribute to the local cosmic ray flux. 
The local supernova rate is about S* = 30 Myr^^ kpc~^. In 
standard cosmic ray diffusion models (e.g. Webber, Lee and 
Gupta 1992) the range of the cosmic rays is linked to the halo 
size ~ 3 kpc and the life time is identical to the escape 
time T ~ 20 7"°-^ Myr. The number of SNR contributing to 
the local cosmic ray (proton) flux is then N = tt t S ~ 
170007^°'^. In other words, while at a proton energy of 50 
GeV some 1600 SNR would contribute, at 100 TeV the num- 
ber would be down to 17 SNR. Even if this is a rather rough 
estimate (using an averaged diffusion coefficient for halo and 
Galactic disc and neglecting reacceleration), it shows clearly 
that the number of SNR contributing to the local CR flux be- 
comes small at energies of some 100 TeV in which case the 
local cosmic ray spectrum could strongly deviate from the av- 
erage spectrum in Eq. ^ A similar effect affects the cosmic ray 
electron spectrum at much smaller energies of 100 GeV, which 
is the basis for the inverse Compton models of the GeV excess 
(Pohl &Esposito |1998D . 

Added to this would be the effect of deviations from power 
law behaviour in the cosmic ray source. We have found that 
there must be a cut-off at proton energies not higher than a few 
hundred TeV. This can be established as sharp cut-offs, possi- 
bly with a distribution of cut-off energies, or also in the form 
of a spectral steepening at somewhat smaller energies, which 
behaviour is predicted in mod els of non-linear shock acceler- 
ation (e.g. Baring et al. 1999). As a result the average spec- 



trum at higher energies would be softer than what is given in 
Eq. ^ which would not change the expected GeV 7-ray spec- 
trum, but the local cosmic ray proton spectrum above 10 TeV, 



which is observed to be slightly softer than that obtained by 
direct measurements at lower energies (Asakimori et al. 1998; 
Amenomori et al. 2000). All in all the lack of dispersion in the 



local cosmic ray proton spectrum does not seem to contradict 
the hypothesis presented in this paper 

Is the high energy limit in the cosmic ray proton source 
spectra compatible with the observed all-particle spectrum and 
the composition near the knee? We have seen that the local 
cosmic ray spectrum near the knee would be produced by few 
nearby SNR which would cause a structured spectrum near 



the knee (Erlykin & Wolfendale 1997a, 



1997b 



1998a, 1998b) 



We have not investigated the composition and in particular 
not studied the possibly different source spectra for different 
species. Therefore we feel not well equipped to discuss the is- 
sue of cosmic ray composition near the knee. The fact that a 
softening in the cosmic ray source spectra is required may in- 
dicate a potential problem with the notion that SNR accelerate 
cosmic rays up to the knee, but a more detailed and careful 
study will be needed to address this issue. 

Since there is observational evidence for a dispersion in the 
spectral indices of the cosmic ray spectra in SNR, the com- 
posite cosmic ray spectrum must be curved, if the cosmic rays 
are predominantly produced in SNR. We have shown that if 
the dispersion is as strong as a = 0.25, its effect on the in- 
terstellar cosmic ray spectrum would explain the GeV excess 
in the diffuse Galactic 7-ray spectrum. If the actual dispersion 
is weaker than this, it would still contribute to the GeV excess 
and therefore should not be neglected. There are other viable 
models for the GeV excess and it may well be that it is a com- 
bined effect of the spectral dispersion and, e.g., a hard inverse 
Compton spectrum (Pohl & Esposito 1998). 



Appendix A: Approximative caicuiation of pion 
spectra 

We consider production of vr^s by cosmic rays via the reac- 
tion pH-p— > tt" + X, where X stands for any other particle, as 
described in Mannheim & Schlickeiser (1994). The total cross 
section for pion production is (Tpp ~ crp,inci — 3 • 10^^^ cm^. 
For charged pions the mean multiplicity has an energy depen- 
dence C ~ 2Ep^^ (GeV) where Ep is in GeV. Due to isospin 
symmetry ^ttO — 5 ^7r±- Laboratory measurements of the mean 
pion energy in p-p and p-a interactions indicate that a constant 
fraction, about 30%, of the incident kinetic energy of protons 
goes to pion energy. This is in agreement with Fermi's theory 
of pion production, in which a thermal equilibrium is assumed 
in the resulting pion cloud and in which the mean Lorentz fac- 
tor of pions in the laboratory frame (where one of the protons 
is initially at rest) is given by 



7^ 



73/4 



(A.l) 



where 7^ denotes the Lorentz factor of the incident proton. 

We approximate the differential cross-section by the exper- 
imentally measured inclusive cross-section multiplied with a 
(5-function centered at the mean pion energy ( A. 1 ), that is 



'^k,pp,pa {l-^ 1 7p) 

k 



^<^pp,pailp) ^il 
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This (5-function approximation and the use of the same mean 
pion energy ( A. 1 ) in p-p and p-a colUsions do not introduce ap- 
preciable error, provided the incident proton is highly relativis- 
tic so that the scaling relations apply. At energies much higher 
than 10"* GeV Eq. (A.2) overestimates the number of pions, but 
at these energies 7-rays are also efficiently produced by kaons, 



so that the use of Eq. (A.2) will not lead to a significant overes- 



timate of the 7-ray source function. (Dermers (1986) formula 



yields an inclusive cross section 20% smaller than that given 



by Eq. (A.2), with an av erage pion energy roughly 50% higher 



than given by Eq. ( |A.2| ), the energy distribution of secondary 
pions going up almost to the kinetic energy of the i nciden t 
proton, overpredicting pions at higher rapidities (Mori 1997 ). 
Using Dermers formula, the cut off in the 7-ray spectrum be- 
comes smoother so the required cut off energy in the CR proton 
spectrum is smaller than in our conservative calculation.) 
For the pion power we then obtain 

P(7^,7) = l.iPc-i^[nHi{r) + 2nHAr)]i^lp{l) 

X ,5(7, - 7^/-*) H[-i ~ ^th] (A.3) 

where the factor 1.30 accounts for the known chemical com- 
position of the interstellar medium, ^th is the hadron threshold 
Lorentz factor The pion source function then is 

Q7r(77r;T") 



d7 A^(7,r)P^n(7^,7) 

1-64 , . X „ 
[nHi{r)+ 2nH2{r)\ x 



/•OC 



d7iV(7,r)/3C<J7)5(7^-7') (A.4) 



where the factor 1.26 accounts for the contribution of a-p, 
a — a collisions and collisions of higher metallicity cosmic 
rays, and is derived using the ratio of the respective inclusive 
cross-section to the p-p inclusive cross-section and the known 
elemental composition of the cosmic rays. 

For a separable differential number density distribution of 
cosmic ray hadrons N{'-f,r) — NQ{r)N{j) at position r = 
(/, b, r) in the Galaxy and a spectrum as given in Eq. (Q) with a 
high energy cut-off at 7max we obtain for relativistic particles 
(/3- 1) 



TTlr, 



^ No{r) [uHiir) + 2nH,{r)] -fl^^ x 



pp ^ 



for 1 « 7. < 7m4x 



(A.5) 



so th at wi th the assumptions made above, the pion source func- 
tion (|A.5l) reduces to 



^ , N 6.6 • 10 ,^ , X r / N „ / M 

Ml^;r) ~ No{r)[nHi{r) + 2nH,{r)] 



X 7,r 



^-^ + ^ ln(7.) 



for l«7-<7n/4 



(A.6) 



77r — 7max 

Eq. ( |A.6[ ) can then be used to calculate the 7-ray source func- 
tion in Eq. (^. 



Acknowledgements. Partial support by the Bundesministerium fiir 
Bildung und Forschung through the DLR, grant 50 OR 0006, is grate- 
fully acknowledged. 



References 

Abramowitz, M., Stegun, I. A., 1972, Handbook of Mathematical 
Functions, National Bureau of Standards, Washington 

Aharonian, F, Akhperjanian, A., Barrio, J., et al., 2001a, A&A 
370,112 

Aharonian, F, Akhperjanian, A., Barrio, J., et al., 2001b, A&A 375, 
1008 

Allen, G. E., Keohane, J. W., Gotthelf, E. V., et al., 1997, ApJ 487, 
L97 

Amenomori, M., Ayabe, S., Caidong, et al., 2000, Phys. Rev. D62, 
112002 

Asakimori, K., Burnett, T. H., Cherry, M. L., et al., 1998, ApJ 502, 
278 

Atoyan, A. M., Aharonian, F A., Tuffs, R. J., & Volk, H. J., 2000, 
A&A 355, 211 

Baring, M. G., Ellison, D. C, Reynolds, S. J., Grenier, I. A., & Goret, 

P, 1999, ApJ 513, 311 
Blandford, R. D., & Eichler, D., 1987, Phys. Rep. 154, 1 
Borkowski, K. J., Rho, J., Reynolds, S. P, & Dyer, K. K., 2001, ApJ 

550, 334 

Brecher, K., & Burbidge, G. R., 1972, ApJ 174, 253 
Buckley, J. H., Akerlof, C. W., Carter-Lewis, D. A., et al., 1998, A&A 
329, 639 

Capella, A., Sukhatme, U., Tan, C.-I., & Tran Thanh Van, J., 1994, 

Phys. Rep. 236, 227 
Clark, D. H., & Caswell, J. L., 1976, MNRAS 174, 267 
Dermer, C. D., 1986, A&A 157, 223 

Engel, R., Ranft, J., & Roesler, S., 1997, Phys. Rev. D 55, 6957 
Erlykin, A. D. & Wolfendale, A. W., 1997, APh 7, 1 
Erlykin, A. D. & Wolfendale, A. W., 1997, APh 7, 203 
Erlykin, A. D. & Wolfendale, A. W., 1998, APh 8, 265 
Erlykin, A. D. & Wolfendale, A. W., 1998, APh 8, 283 
Esposito, J. A., Bertsch, D. L., Chen, A. W., et al., 1999, ApJS 123, 
203 

Ferrari, A., Sala, P R., Ranft, J., & Roesler, S., 1996a, Z. Phys. C 70, 
413 

Ferrari, A., Sala, P R., Ranft, J., & Roesler, S., 1996b, Z. Phys. C 71, 
75 

Green, D. A., 2001, in High energy gamma-ray astronomy, eds. F. 
A. Aharonian and H. J. Volk, AIP Conference Proceeding 558, 
59,2001 

Hartman, R. C, Bertsch, D. L., Bloom, S. D., et al., 1999, ApJS 123, 
79 

Hunter, S. D., Bertsch, D. L., Catelh, J. R., et al., 1997, ApJ 481, 205 
Koyama, K., Kinugasa, K., Matsuzaki, K., et al., 1997, PASJ 49, L7 
Koyama, K., Petre, R., Gotthelf, E. V, et al., 1995, Nature 378, 255 
LeBohec, S., Bond, I. H., Bradbury, S. M., et al., 2000, ApJ 539, 209 
Mannheim, K., & Schlickeiser, R., 1994, A&A 286, 983 
Milne, D. K., 1979, Austr. J. Phys. 32, 83 

Moffett, D. A., Goss, W. M., & Reynolds, S. R, 1993, AJ 106, 1566 

Mohring, H.-J., & Ranft, J., 1991, Z. Phys. C 52, 643 

Mohring, H.-J., Ranft, J., Merino, C, & Pajares, C, 1993, Phys. Rev. 

D47, 4142 
Mori, M., 1997, ApJ 478, 225 

Muraishi, H., Tanimori, T, Yanagita, S. et al., 2000, A&A 354, L57 
Pohl, M., & Esposito, J., 1998, ApJ 507, 327 
Ranft, J., & Roesler, S., 1994, Z. Phys. C 62, 329 



Busching, Pohl, Schlickeiser: Excess GeV radiation 



7 



Ranft, J., Capella, A., & Tran Thanh Van, J., 1994, Phys. Lett. B 320, 
346 

Reynolds, S. P., & Keohane, J. W., 1999, ApJ 525, 368 
Roesler, S., Engel, R., & Ranft, J., 1998, Phys. Rev. D 57, 2889 
Slane, P., Gaensler, B. M., Dame, T. M., et al., 1999, ApJ 525, 357 
Sreekumar, P., Bertsch, D. L., Dingus, B. L., et al., 1998, ApJ 494, 523 

Stecker, F. W., 1973, ApJ 185, 499 

Tanimori, T., Hayami, Y, Kamei, S., et al., 1998, ApJ 497, L25 
Webber, W. R., Lee, M. A., & Gupta, M., 1992, ApJ 390, 96 



